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I,  Introduction  and  ■uamarv. 

Among  tho  bnilc  concoptt  In  tho  thoory  of  OADpllng  Inspoctlon  cho  producer'* 
and  the  consumer'*  risk*  are  the  most  widely  used  for  charscterltln3  sysceia. 
of  sampling  plans.  It  seems  therefore  strange  that  a  comprehensive  theory  based 
on  these  concepts  does  not  exist.  The  purpose  of  the  present  paper  is  to  present 
such  a  theory  for  the  ease  of  single  sampling  by  attributes.  The  theory  naturally 
covers  some  well-known  results,  but  old  as  well  as  new  results  are  derived 
by  a  coisBon  method  and  compared  within  the  same  model.  The  requirements  defining 
a  system  of  sampling  plans  are  usually  of  such  a  nature  that  no  explicit  solution 
exists  for  the  sample  sise  and  the  acceptance  number.  We  shall  therefore  supple¬ 
ment  the  exact  (implicit)  solutions  by  asymptotic  solutions  which  give  a  better 
insight  into  the  basic  properties  of  the  systems. 

Let  there  be  given  two  quality  levels,  p^  and  P2,  ^  P2*  ^  sampling  plan, 

(n,c),  n  denoting  the  sample  eiae  and  c  the  acceptance  number,  the  operating 
characteristic  is  defined  as  P(p)  ■  B(c,n,p),  where  B(c,n,p)  denotes  the 
cumulative  binomial  distribution.  The  producer's  risk  is  then  Q(pj^)  •  l-P(Pj^), 
and  the  consimier's  risk  is  P(p2)  •  These  risks  give  the  probabilities  of  %rrong 
decisions  under  the  assumption  that  p^^  represents  acceptable  and  P2  rejectable 
quality. 

We  shall  furthermore  aesiaae  that  the  consequences  of  wrong  decisions  are 
comnsnsurable  and  measurable  and  that  the  average  "loss"  from  using  a  given 
sampling  plan  may  be  expressed  as  a  linear  combination  of  the  two  risks, 

7jQ(Pi)  72^(P2)  **7'  *  Bayesian  point  of  view  7^  equals  the  product 

of  the  prior  probability  of  p^  and  the  corresponding  decision  loss. 

If  the  costs  of  sampling  inspection  are  proportional  to  the  sample  size  and 
we  sample  a  lot  of  size  N  we  may  therefore  write  the  average  costs  in  the 
standard  form 

R(N,n,c)  -  n  +  (N-n)(7jQ(pp  +  72PCP2)) 

where  7^  1  0  and  72  k  0.  We  shall  use  this  cost  function  in  the  comparison  of 
the  various  systems  of  sampling  plans. 

In  section  2  it  is  shown  how  B  may  be  interpreted  as  the  average  costs  when  (1) 
samples  are  drawn  without  replacement  (2)  from  lots  produced  under  binomial 
control  but  with  a  process  average  varying  at  random  from  p^^  to  P2,  i.e.  the 
prior  distribution  of  lot  quality  is  a  double  binomial,  and  (3)  costs  are 
linear  in  the  number  of  defectives,  the  sample  sise  and  the  lot  size.  Apart 
from  a  term  of  order  *"  the  function  R  may  thus  be  interpreted  as  the  usual 
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"risk"  in  decision  theory.  The  model  covers  rectifying  as  well  as  non- 
rectifying  inspection. 

In  sections 3  and  A  asymptotic  expansions  are  derived  for  the  producer's  and 
the  cooeumr's  risk  and  for  their  ratio. 

Section  5  gives  a  minimization  theorem  from  which  the  relation  between  lot 
site  and  sample  site  may  be  found. 

In  the  remaining  part  of  the  paper  we  discuss  ten  systems  of  sampling  plans 
defined  as  follows: 

(1).  Bayesian  plans,  i.e.  plans  mlntmizlng  R. 

Restricted  Bayesian  plans,  i.e.  plans  minimizing  R  under  some  suitably  chosen 
restriction  on  the  operating  characteristic,  viz. 


(2). 

Min  R  for  Q(Pj)  ■  a  or  P(p.,)  ■  P. 

\  V2/ 

(3). 

Min  R  for  P(Pq)  •  1/2  where  p^-  ^log^ )  y 

(A). 

Min  R  for  Q(p^)  ■  a/N  or  P(P2)  "  P/N. 

(5). 

Min  R  for  P(P2)/Q(pp  "  P. 

Plans  defined  by  two  risks,  vis. 


(6). 

Q(Pj)  =*  a/N  and  Pfp^) 

-  P/N. 

(7). 

Q(P^)  "  a  and  PiPj)  “ 

P/N  (or  P(P2)  ■  P  and  Q(pp  -  a/N) 

(8). 

P(Pq)  -  1/2  and  Q(p^) 

-  u/N  (or  P(P2)  -  P/N) . 

(9). 

Q(Pj^)  ^  a  and  P(p2)  =* 

P. 

Finally  we  consider  percencap.e  InspccLlon  defined  as 
(10)  .  n  =»  uN  and  c  ■  p^n. 

In  all  these  definitions  and  p  represent  suitably  chosen  positive 

constants  which  may  be  different  from  case  co  case. 

For  each  system  of  sampling  plans  it  is  shown  how  the  exact  solution  nay  be 
obtained  and,  since  this  solution  is  an  implicit  one,  an  explicit  solution 
is  given  as  an  asymptotic  expansion  for  N  ">  oo. 

The  advantages  of  the  asymptotic  solution  are  that  (1)  it  clearly  shows  how 
the  sampling  plan,  the  two  risks,  and  the  costs  depend  on  the  parameters,  (?.) 
it  gives  good  approximations  to  the  exact  solution  even  for  quite  small 
values  of  c  (normally  sufficiently  accurate  for  c  2  2),  (3)  it  may  be  used 
for  developing  interpolation  and  extrapolation  formulas  in  connection  with 
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"master  cables"  of  the  exact  solution,  and  (A)  ic  shows  the  sensitivity  of  the 
solution  with  respect  to  changes  of  the  parameters.  We  shall  here  mainly 
discuss  the  first  of  these  points.  With  respect  to  the  other  three  and  a  de¬ 
tailed  discussion  of  tables  the  reader  is  referred  to  [9],  [11],  and  [12]. 


The  systems  defined  by  (1)  -  (8)  fall  into  two  classes  depending  on  whether 
both  risks  are  0(H  or  one  of  the  risks  Is  constant  and  the  other  is  0(N  . 

The  first  class  -tontaLns  systems  (1),  (3)^  (4),  (5),  (6),  and  (8).  It  is  proved 

that,  asymptotically  the  relation  between  acceptance  number  and  sample  size 

•1*2  1  *1  *2 
has  the  form  c  ■  p^n  -t**  a^n  +  0(n  ),  that  In  N  *  (p^n  +  jln  n  +  Kj^+  K2n  +0(n 

which  by  inversion  determines  n  as  a  function  of  N,  and  that  R  ■  n  +  h^n  Vo(n  ^), 

the  constants  p^  and  9^  being  the  same  in  all  canes  (depending  on only), 

whereas  the  remaining  constants  are  found  as  functions  of  the  parameters  in  the 

model  and  the  restriction.  This  means  that  n  »  0(ln  N)  and  that  the  average 

decision  loss,  (N-n)  (yj^Q(Pj^)  +  tends  to  a  constant  (because  Q(pp  and 

^(p^)  are  0(N"S). 


The  second  class  consists  of  systems  (2)  and  (7).  Ic  is  proved  that  asymptotical¬ 
ly  the  relation  between  acceptance  number  and  saraplc  size  has  the  form 
^  1  — ''/2  —3/2 

c  *■  p  n  +  Z  a,n  ~  +  0(n  ^),  p  representing  the  quality  level  having 

^  i.  ■  I  ^ 

’  ^  -3/2.  .  .  . 

+  0(n  ),  whLch 


1/2, 


a  constant  risk,  that  In  N  »  (p.n  +  ^In  n  +  Z  K.n 

J  ^  i=l 

determines  n  as  a  function  of  N,  and  that  R  --  6N  +  (l-6)n  +  6^+  0(n 
the  constants  9^  and  5  being  the  same  in  the  two  cases.  Because  of  the  constant 
risk  all  relations  are  considerably  more  complicated  than  for  the  first  class 
and  R  becomes  0(N)  instead  of  0(ln  N) .  For  large  lots  it  must  therefore  be 
seriously  considered  whcchcr  it  is  reasonable  to  use  a  system  with  a  fixed 
consiimer^s  or  producer's  risk  and  correspondingly  high  costs  as  compared  with 
a  system  having  decreasing  risks. 


The  system  with  both  risks  fixed  and  the  system  with  percentage  inspection  both 
icad  to  R  ■  0(N)  and  asymptotically  they  have  the  same  costs  for  p  =  7^a  + 

The  system  with  fixed  risks  uses  a  fixed  sample  size  so  that  the  decision  loss 
becomes  of  order  N  whereas  percentage  inspection  has  n  ■  0(N)  and  a  decision 
loss  of  order  e 


For  systems  (2)  and  (3)  we  have  the  Important  result  that  asymptotically  n 
depends  on  (Pj,P2)  and  NX  only  where  X  Is  a  function  of  7^t72f  parameter 

In  the  restriction.  It  cherefore  suffices  to  tabulate  n  as  a  function  of  N  for 
X  ^  1,  say,  and  use  this  table  for  X  4  1  with  N*®  NX  as  argument. 

Writing  B(c,np)  for  the  '.umulative  Poisson  distribution  corresponding  to  B(c,n,p) 
we  have  under  Poisson  conditions  exactly  and  under  binomial  conditions  for  small 


.  4  - 


(P11P2)  approximately  that 

Rpj^  «=  m  +  (M-fa)(7^(l-B(c,m))  +  72*(c,na)) 

%rhore  a  ■  np^^  M  »  Np^,  and  r  "  analogous  to  those  stated  above 

are  therefore  valid  in  terms  of  M,n,  and  c  for  given  (r, 7^,72)^ 
one  parameter.  It  follows  that  under  binomial  conditions  %re  have  for  systems 
(1)>  (2),  (3),  and  (5)  the  foll(nring  "proportionality  law":  The  sarapling  plan 
corresponding  to  (N;Xp^,Xp2)  is  approximately  equal  to  (n/X,c)  ifhcrc  (n,c)  is 
the  plan  corresponding  wO  (NX,PpP2). 

For  systems  (6)  ■■  (9)  we  have  similarly:  The  sampling  plan  corresponding  co 
(N^Xp^^Xpj)  is  approximately  equal  to  (n/X,c)  where  (n,c)  is  the  plan  correspond¬ 
ing  to  (N^p^iPj). 

These  theorems  greatly  reduce  the  tables  necessary  fur  applications  of  the 
systems. 

It  is  sho%m  that  the  restricted  Bayesian  plans  with  both  risks  decreasing  and 
the  corresponding  plans  based  on  two  risks  all  have  an  economic  efficiency  tend¬ 
ing  to  1  for  M  00  as  compared  to  the  Bayesian  plans.  (The  efficiency  of  plans 
having  at  least  one  risk  fixed  tends  to  zero) .  This  result  means  that  wrong 
values  of  the  weights  of  the  prior  distribution  and  wrong  values  of  the  cost 
parameters  have  a  secondary  influence  on  the  efficiency  which  tends  to  1  if 

only  (p.,p«)  are  correct.  If  also 

^  #  # 

wrong  values  of  (PpP2)>  (P1/P2)  say,  are  used  fer  finding  the  plans,  then  the 

“  #  # 

efficiency  tends  to  e,  0  <  e  <  1,  if  and  only  if  ^  Pj^  ^  p2  ^  P2^  otherwise 
the  efficiency  tends  to  0. 

The  present  model  leads  to  a  constant  ratio  of  the  producer's  and  the  consumer’s 
risk  for  the  (Bayesian)  sampling  plans.  This  provides  a  motivation  for  the  rule 
of  thumb  suggested  by  Lehmann  [14]  for  obtaining  a  reasonable  balance  between 
the  probabilities  of  errors  of  the  first  and  second  kind  in  testing  the  hypothesis 
p  “  Pj  against  the  aiternative  p  ■  P2* 

2.  The  model. 

Let  N  and  n  denote  lot  size  and  sample  size  and  let  X  and  x  denote  number  of 
defectives  in  the  lo.  and  the  sample,  respc  ively.  The  acceptance  number  is 
denoted  by  c 

Consider  the  following  linear  cost  function 

nSj^+  XS2+  (N-n)A^+  (X-x)A2  for  x  S  c 

h(X,x,N,r.,c) 

I 

inS^+  XS2+  (N-n)R  +  (X-x)R2 


for  X  >  c 


(1) 


-  5 


r.nd  let  tho  (prior)  distribution  of  lot  quality  be  f^^CX).  The  average  costa  then 
become 

K(N,n,c)  E  E  h(X,x)  p(X,x) 

X  X 


where 


p{X,x)  . 


8«  (x)  giving  the  (marginal)  probability  of  getting  x  defectives  in  the  sample, 

N,n 

i.c. 


N-n 


'  yr  O 


(2) 


This  compound  hypergeomctrlc  distribution  has  been  discussed  in  [s]  %dierc  Che 
following  results  were  proved: 


and 


ri 

x<') 

n 

.n<^) 

eI 

(X-x) 

(N-n) 

c 

1  s 

Z 

x-o 

K(N 

^n,c) 

,(3) 


»  E 


ry(r+8)' 


i 

'-N 


(r) 


(r) 


I 


(n+r) 


(r) 


(X) 


’N,n 


'N,n'’'  “  ®N,n 

K(N,n,c)  -  nS^  +  E{x)S2  +  (N-n) 


(3) 


(A) 


n 


+  /‘2  ^  j^(x)E{X-x|x)  +  E  g„  ^(x)E{X-x|x].  (5) 

x«=o  *  x-c+l  ' 


The  last  two  terms  of  (5)  may  be  simplified  by  using  (''O  for  r  «  1  which  gives 

The  average  costs  are  obviously  a  rather  complicated  function  of  n  and  c  and 
in  general  also  of  N.  It  is  easy  to  see,  however,  that  the  necessary  and  sufficient 

condition  for  K  to  be  a  linear  function  of  N  for  all  values  of  (n,c)  and  all 

values  of  the  cost  parameters  is  that  g  (x)  is  independent  of  N. 

Cl  ^  n 

From  g^s  independence  of  N  follows  that  E(x^*^Vn^^^)  is  independent  of  N  so 
that  (3)  gives 


ECX^’^^)  -  a 


for  N  J  r. 


(3) 


say,  where  does  not  depend  on  N, 
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To  find  the  class  of  distributions  satisfying  (G)  wc  introduce  the 

limiting  cumulativci  distribution  of  X/N  defined  by 


W 

H(p)  -  Urn  I  £  (X)  . 
>  OD  X"0 


(7) 


Proceeding  as  in  [o]  we  find  that 

1 

a  »  /  p’^dlKp) 

0 


and 


f^(X) 


-  (”) 


I 

/ 

0 


X  N-X 

p  q 


dW(p) 


(G) 


It  follows  that  g„  (x) 

Nf  n 


f^(x),  i.e.  the  prior  distribution  is  "reproduced"  by 


hypergeomctric  sampling.  discussion  of  reproducible  distributions  has  been  given 
in  [o]). 


We  shall  therefore  get  particularly  simple  results  by  .limit ir^_thc_2r^r 
distributiona  to  the  class  of  distributions  given  by  (S).  called  mixed  binomial 
distributions. 

Writing  b(x^n,p)  «  *  and 

c 

P(p)  -  B<c,n,p)  «  E  b(x,n,p)  (9) 

x*o 


we  get 


and  from  (5) 


where 


x-fl 


^  b(xfl,n+l,p)  *  p  b(x,n,p) 


1 

K(N,n,c)  -  /  K(N,n,c,p)dW(p) 

0 


K(N,n,c,p)  -  n(Sj+S2P)  +  (N-n) ^^+A2P)P(p)  +  (Rj+  R2P)Q(P) 


(10) 


.  (11) 


The  assumption  of  a  mixed  binomial  prior  distribution  means  that  each  lot  is 
produced  under  binomial  control  and  that  the  process  average  varies  at  random 
from  lot  to  lot  according  to  the  cumulative  distribution  function  W(p). 
Correspondingly  the  average  costs  (11)  represent  on  average  over  all  lots  with 
a  given  process  average,  i.e.  a  conditional  average,  and  (10)  gives  the  over-all 
average. 


Besides  giving  the  exact  average  costs  for  a  mixed  binomial  prior  distribution 
(10)  and  (11)  may  be  interpreted  as  giving  an  approximation  to  the  average  costs 
for  laigc  N  for  any  prior  distribution  satisfying  (7).  This  follows  from  the 
facts  that  the  hypergeomctric  distribution  tends  to  the  binomial  for  N  — >  oo, 
n  — 00,  n/N  — >  0  and  X  *  Np,  p  fixed,  and  that  x  np  +  0(i/h),  so  that  (1) 


I 

*1 


becomes 


r  n(Sj  +  $2?)  +  (N  -  n)(Aj  +  A2P) 

t  n(Sj  +  S2P)  +  (K  -  n)(Rj  +  R2P) 

disregarding  terms  of  order  vn.  The  limit  theorems  derived  in  the  following 
on  the  basis  of  (10)  and  (11)  are  therefore  valid  in  general. 

The  model  may  easily  be  generalized  from  a  linear  to  a  polynomial  cost  function. 
Consider  for  example  the  term  (X  -  x)A2  %fhich  gives  the  contribution 

1 

(N  -  n)A«  /  p  P(p)  dW(p) 

0 

to  the  average  costs.  Introducing  instead 


m 


Z  A2^  (X  -  X) 
v*=l 


(v) 


we  find  by  using  (A)  the  following  average  costs 

v-1  ‘^N,n+v 


(x+v). 


The  condition  for  the  average  costs  to  be  a  polynomial  in  N  is,  as  above,  that 
g  docs  not  depend  on  N.  For  a  mixed  binomial  prior  distribution  we  get 

ro  /  \  1 

S  A  (N  -  n)'"-'  /  p^'pfp)  dW(p). 
v-1  0 

Treating  all  six  terms  of  (1)  analogously  we  find  the  generalized  average  costs 
are  given  by  (10)  if  we  replace  (11)  by 


m 


K(N,n,c,p)  -  I  +  S^y) 

v*l 


m 


+  I  (N-n) 

V"1 


<v) 


(Aiv'^2vP'^)P(p)  ■*' 


Another  generalization  which  is  easily  carried  out  consists  of  replacing  (8)  by 


f^(X) 


I  X  N-X 

L)  /  p  q 
^  0 


<IWj^(p) 


(13) 


where  the  comulativc  distribution  function  W  (p)  dep'tnds  on  N.  This  will  only 

w 

result  in  a  corresponding  change  of  (10). 


In  the  following  we  shall,  however,  mainly  discuss  the  cost  function  defined 
by  (10)  and  (11). 


To  siiapllfy  the  no:aclon  we  Introduce  the  three  cost  functions 


ka(p)  -  A^+  AjP,  ky(p)  -  R2p>  kg(p)  “  8^+  S^p,  (lA) 

and  che  corresponding  averages  k^/  defined  by 

1 

k  ■  /k(p)dW(p)  .  (15) 

0 

If  the  equation  k  (p)  »  k  (p)  has  a  solution  p  "  p  In  the  interval  (0^1)  we 

a  r  r 

define 

**r  1 

k  -  /  k  (p)Sw(p)  +  /  k  (p)dW(p)  (16) 

0  ^  p/ 

%#hlch  represents  the  average  costs  per  iteu  when  all  lots  from  processes  with 
p  $  p^  are  accepted  and  all  other  Iocs  arc  rejected.  The  fraction  defective 
p^  is  called  the  (economic)  break-even  quality. 

Defining  the  standardized  form  of  (10)  as 


R(N,n,c)  -  (K(N,n,c)  -  Nk  )/(k  -  k  )  (17) 

m  s  m 

we  find 

A--R  ^r  I 

R  -  n  +  (N-n)  i  /  (p^-p)Q(p)dW(p)  +  /  (p-p^)P(p)dW(p)  ,  (18) 

'am-  p  ^ 

Using  k^-  k^  as  "economic  unit"  the  two  terms  of  (18)  represent  the  costs  of 
sampling  inspection  and  the  average  decision  losses^  respeccively . 

In  the  remainder  of  the  paper  ic  will  be  assumed  that  the  prior  distribution 
is  a  double  binomial  distribution,  or  as  a  limiting  case,  a  single  binomial. 

The  double  binomial  di.sCribution  is  a  weighted  average  of  two  binomials  with 
parameters  p^  and  Pj^<  P2,  and  weights  w^  and  W2,  w^-f  the  process 

average  has  a  two-point  distribution. 


The  standardized  average  costs  may  then  be  written  as 


R(N,n,c)  -  n  +  (N-n)(r^Q(pp  +  72P(P2)) 


(19) 


where 


7,“  ''i(k_(p,)-k  (p,))/(k  -  k  )  and  7  -  w  (k  (p,) "k  (p  ))/(k  -k  )  .  (20) 

I  trial  SQ  c  i  a  i  t  2  s  n. 


It  will  be  noted  that  the  average  decision  loss  per  item  is  a  linear  combination 
of  the  producer's  risk.  Q(pj^),  and  the  consumer's  risk.  P(P2) 

The  standardized  costs  of  accepting  or  rejecting  all  lots  without  inspection 
are  R^*  N72  and  R^=  N7^,  respectively. 
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3.  An  asymptotic  cxpanoicn  for  the  binomial  distribution. 

There  exist  several  asymptotic  expansions  for  1-B(c,n,p),  c  ■  [np^]  and  p  <  p^/ 
for  n  ~>  see  for  example  Blackwell  and  Hodges  [1]  and  Brockvell  [2].  We  need^ 
however,  an  expansion  under  the  assumption  chat  c/n  ■  pQ*f  c,  c  — >  0  for  n  — >  «, 
and  shall  use  the  same  method  as  Brockvell  to  prove 


Theorem  1 .  Le  t 


-i/2  -5/2 

c/n  ”  ^  d-  0(n  ) . 

”  1-1  ^ 


(21) 


For  p  >  Pq  we  have 


B(c,n, p) 


exp  ■|^-n((*»(pQ,p)  +  Z  h^(pQ,p)n  ^^^40(n  (22) 


where 


|P'Pol'2rtnpQqQ 


PQ 


b  *  a. In  ^  , 

1  ‘  <loP 


''o 

1%.?)  ‘  Pq'-'T'^  ''o‘''T  ' 

2 


Pq** 

b«-  a«,ln — -  + 


2  2  loP  ^PoS) 


Pfli  ^  ^<VPo> 

b-*’  a*  In -  -  - +  ^  _  -  - r  , 

^  ^  ^P  P-Po  ^^Po-0 


b, “  a,  In-^  - 

lo 


a  %A  a  *  •  • 

^  ^  P"p- 


0  2(p-Pq) 


^  ^-PqIq  <’lo•Po>°^‘*^P2>  ,  <Po*^>''t 
‘^PO-^  MPgqo)^  I2(Po<lo>^ 

Foi  p  <  pQ  the  same  expression  is  valid  for  1-B(c,n,p). 


Proof.  Writing  c  »  nh  and  using  Stirling's  formula  we  get 


ln(") 

so  that 
In  b(c^n,p) 
where 


i  ln(2nnh(l-h))  -  n(h  In  h  +  (1-h) In(l-h))  -  +  0(n’^) 


'-|ln(2ffnh(l-h))  -  n(p(h,p)  -  '}~Jnh(i-h) 
c;;(h  i)  -  h  In—  +  (l-h)ln-^^  . 


(23) 


Expanding  (p(h,p)  in  a  Taylor  series  around  p  and  inserting  the  given  expression 


for  h-p^  we  find 


(p(h,p)  -  (Pq+  a^cp^n  ^®2'^l‘*'  2®  1^2^”  ^“3‘^l‘*"  *1*2‘^2'*'  6®l‘^3^"'^^^ 

+  (a^cpj  +  (|n^+  ajaj)(p2+  ^  +  0(n’^^^ 

where  cp^-  (d^(p/ Jh^) .  Similarly  wc  have 
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2.  2.  2 


l-h(l-h)  ^  -3/2, 

I2nh(l-h)  "  12npjjqQ  ^  ' 

Combining  these  expressions  we  get  an  expansion  for  In  b(c,nyp)  wl  eh  remainder 
eenn  of  order  n”  '  . 

Expressing  the  binomial  by  the  incomplete  DcLafunction  we  have 

B(c,n,p)  ■  (n-c)(")  /(1-x) ‘^x"'‘^‘^dx. 

0 

Changing  Che  variable  the  integral  becomes 

.  n 

00  00  /  n  \  00 

/  <y-i) /  o 

1/,  1/,  ^  y  '  y 

where  s  ■  -h  In(y-l)  -!•  In  y  and  ffz)  *  y  ^dy/dz  Integrating  by  parts  leads  to 
00  ^nz 

/  e"''*f(z)dz  -  K  °  (f(z  )  +  ^f^(2  )  +  0(n’^). 

^  n  0  n  0 

As  Zq  is  found  from  z  for  y  -  1/q  we  have  z  -h  In  p  -  (l-h)ln  q  and 

-nz^ 

Oj  c  n-c. 
c  /n  •  p  q  /n.  From 


we  get 


f(z)  .  (y-l)/(y(l-h)-l) 


f'(z)  -  -hy(y-l)/(y(l-h)-l) 


so  that 


n(p-h)  ^ 


Combining  the  results  obtained  we  have 

B(c,n,p)  =  b(c,n,p)^^'^^P  Tl - — r  +  0(n’^) 

P"^  ^  n(p-h)^ 

Expanding  the  logarithm  of  the  factor  to  b(c,n,p)  around  p^  we  find 

1„  B(c,n,p)  .  V.  ^  ±'' 

b(c,n,p)  P-Pq  t,n  Vp-Pq  q^^ 

,  i  f  ,  vfibL  .  !2 . 4  ^ 
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which  result  together  with  the  expansion  for  In  b(c,n, p)  after  some  reductions 
lead  to  the  theorem. 

A  similar  procedure  leads  to  the  result  for  l*B(c,n,p). 

It  Is  obvious  that  the  expansion  may  be  continued  by  using  the  sarae  method 
and  that  the  following  terms  are  of  the  same  type  as  those  given  in  theorem  I. 


4.  The  ratio  of  the  consumer's  to  the  producer's  risk. 

From  theorem  1  it  follows  that  both  the  consumer's  risk,  P(P2)  ■  and 

the  producer's  rick,  Q(pp  *  l"B(c,n,Pj^) ,  tend  exponentially  to  zero  for  n  — >  * 
and  p^<  p^<  p^  sin-c  cp(pQ,p^)  >  0  and  epfp^,  p^)  >  0. 

To  discuss  the  ratio  wc  Introduce 

P2<»i  Qi 

6(p)  9(P>P,)-v(PiPo)  =  P  In  -  In  —  , 

1  2  p^q^  q2 

6o=*  ^(Pq)  nnd  6*  =  In(p2q^/p^q2)  •  I^  follows  from  theorem  1  that 

«i(pj-p,)  ^ 


(25) 


In 


P(P2) 

Q(iV 


In  ■  ■  ..  +  no 


/Pj-Pq) 


iQ+Jm>j6'+0j(f+^  ^36'+ 


(Pj-PoXPo’Pi)  / 


+  -  (a,S’+  ;■  - r  +  — 4 - - * - Si-4  -KlCn  " .  (26) 

V  '  (Pj-PqXPo-P,)  2(p  .p„)^  2(p„-p,)^'^ 


sS2(P2-Pi) 


-3/2, 


2  *^0' 


rQ  ... 


n6. 


n 

For  6^4  0  find  that  P(P2)/Q(Pj)  ■  0(e  ),  so  that  one  of  the  risks  tends 

exponentially  faster  to  zero  than  the  other. 

For  6q*  0  and  =  0  wc  have  that  P(P2)/Q(pp  tends  to  a  constant. 

By  means  of  the  above  expression  wc  shall  prove 

/  ,  /  P2*>  A 

Theorem  2.  Let  c  =  nD_+  a»  and  p_»  In  —  ;  In  —  ^ 

-  2  ^0  V  q2//\ 


.  Then 


^  .  P;<Po~Pl>  W  j  ^  _i 
Q(Pl)  “  P](P2"V  ^ 


^-^2  ^  ^'^^2  *^0  *b’^l'  (P2"Po^^ 


+ +  0(n‘^''^)  .  (2 7) 


For  sraall  p^  and  p^  and  02“  -2/3  we  have  approximately  P(Pq)  •  1/2  and  P(P2)  ■  W(Pj)  • 

Proof.  The  first  jiart  of  the  theorem  foll(3W8  from  (26)  for  6^*  a^*  a^"  a^«  0. 

The  second  part  is  found  by  letting  p^->  0  for  fixed  r  =  P2/P^-  Noting  chat 
p^/p^  ->  (r-l)/ln  r  and  introducing 

^1  “  r-l-ln  r  ^2"  r  In  r-r+l  , 

I 

I 
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we  find 


P(P2) 

^PjT 


82+1 


82 


which  for  r  <  20  Is  approximately  equal  to 


P(P2) 

QCPj) 


82+2/3 


The  last  result  rests  on  the  (numerical)  fact  that  r 
for  r  <  20  as  will  be  seen  from  the  following  table. 


1/3 


-3/2, 


1  +  (81+  82)(a2+  Si"  82)  f  +  0(n  ) 


82/81  a  1  and 


(2C) 

82  ^  ^/^ 


Table  of  g^,  g^,  and  r 


1/3 


82/Ci 


r 

81 

82 

81*82 

81+82 

1/3  , 

r  82/81 

1 

00 

00 

0.37 

00 

1.00 

2 

3.23 

2.59 

0.67 

5.C5 

1.00 

3 

2.22 

1.54 

0.  3C 

3.76 

1.00 

5 

1.67 

0.99 

0.  38 

2,3  3 

1.01 

7 

1 .  AG 

0.79 

0.^9 

2.27 

1.02 

10 

1.3A 

0.6A 

0.70 

1.90 

1.03 

15 

1.2A 

C.53 

0.71 

1.77 

1.05 

20 

1.1‘. 

0,46 

0.  73 

1...5 

l.Oo 

Since  the 

equation 

B(c,n,pp)  = 

1/2  Has 

the  solution  c  ■  np 

-  (2  -  Pq)/3  +  g(n’S, 

sec  section  S,  the 

last  part  of 

the  theorem  is  proved. 

Theorem  2 

implies  that  a  system 

of  sampling  plans 

def ined 

(partially)  by  the  relation 

c  .  nPo  + 

will  have  a  ratio 

decreasingly  or  increasingly  to 

a  constant 

according  as  02  is 

larger 

or  smaller 

than-2/3 

respectively.  If  p^  differs 

from  the  value  defined  in  theorem  2  then  PXP2)/Q(Pj)  will  tend  exponentially  to 
zero  or  infinity. . 

From  (26)  we  also  find 

^  ^  y  0 

Theorem  3.  Let  c/n  =  ^  a.n  ^  +  0(n  ).  Then  P(p^)/Q(p,)  is  constant 

i-1  ^  ^  ‘ 

(to  the  order  of  approximation  here  considered),  i.e. 


P(P2) 

QfPj) 


(!2h\  V, 

Pl(P2-Po)  vPiq2^  - 


(2S) 


if  and  only  if  p,  -  (In  — )/ln  -  ,  a,  •*  a_  -  0,  anu 

0  <12  vPi‘l2/  '13^ 


13  - 


a 

4 


l  -Y  °2 

“  Vpo’Pi 


Po-P 


_Poli 


2  ‘'0  (Pq-p^)' 


Pn<l 


0'*2 


(Pj-Pq)' 


(30) 


For  anall  and  p^  we  have  approximately 

c  a  np^  +  a^  -  (3^+  82X02+  2/3)/(npQln  r) 

*2+  2/3 

and  P(P2)/Q(Pj)  a  r  where  r  ■  P2/Pj. 

Prcof,  The  main  result  follows  directly  from  (20)  and  Che  approximation  is  found 
by  similar  considerations  as  under  theorem  2. 


5.  A  minimization  theorem. 

For  use  in  later  sections  we  need 
Theorem  4.  The  minimum  of 


-1/2 

R(n)  "  n  +  (N  -  n)Xn  f(n) 

^  -i/2 

where  f(n)  ■  exp{-n  E  a  n  ),  \  >  0  and  a  >  0,  with  respect  to  n  for  N 

i-0  ^  ^ 

is  equal  to  2 

1  /  “1  “1  "  ^*^0^  -3/2 

min  P  «  n  +  •—  f  1  -  A'  +  ^  *  +  0(n 


a. 


2a^  V^n  /  2 
0  ‘ 


n  being  determined  ns  a  function  of  N  from  the  equation 

1  ^  -l/'’  -3/2 

In(N-n)  «  -  In  \+  -r  In  n  +  n  E  a  n  *"  +  0(n  ) 

i=0 

where  -  -Tj  -  In  a^,  and 

^40  “  "4  -  ‘/'“O  +  • 


(31) 
->  00 

(32) 


(33) 


Proof.  From 


RXn)  »  1  -  xn'^^^f(n)  +  (N  -  n)(Xn"^^^f Xn)  -  jXn'^^^f(n))  -  0 


we  get 


(N  -  n)xn  ^'^^f(n) 


Xn  f(n)  . 


It  follows  that  n  — >  oo  for  N  — 00, 


Since 


.  IlM  +  i- 

f(n)  2n 


1  1  -3/2 

“0  2^  2^;  > 


(34) 


and  1 


\n"^^^f<n) 


1  +  0(e”*^)  we  find 


(N  -  n)Xn*^^^f(n) 


( 


°1  I  -3/2 

2a-  ^  2 


0  \ 


4Con  / 


+  0(n‘^^^) 


which  innediaCely  gives  min  H.  Inking  logarithms  on  both  aides  of  the  last  equation 
we  obtain  the  equation  for  determining  n. 


Inversion  of  this  equation  gives 


*=  X  -  ^In  X  +  02  +  ^3^*  ^^^In  x  +  3^*  ^  ^  ®(*  (^5) 

-1/2  12  3 

where  x  =  In  N,  ,  p2  =  -jp^-  ^2  2'^^  '0  ■**  ' 

13  1  1/2  3  1 

*^3  ~  ^1  ’  0*^1  2^1°2  '  ^’g  ^3  '  4^1  “^0  *  2^1  ^  ^ 

13  2  1  3  1 

3/  “  T  "  “3i  +  "  7!*^  x» 

4  z  ui  Z  Z  04  4  0  Z 


which  determines  the  value  of  n  in  min  R, 

The  importance  of  theorem  4  is  due  to  the.fact  that  the  asymptotic  form  of 

4  X  *  /  2 

R(N,n,c)  after  ellaiaation  of  c  -  np^  +  I  an  is  equal  to  (31),  The  result 

O  1  i 

(32)  then  tells  that  the  minimum  standardized  costs  asymptotically  consist  of 
sampling  inspection  costs  plus  a  term  tending  to  a  constant  which  represents  the 
limiting  average  decision  losses  for  the  remainder  of  the  lot. 


An  important  corollary  is  found  by  noting  from  (33)  that  asymptotically  n  depends 
on  13L  on  ln(NX.),  i,e,  on  the  product  of  the  lot  size  and  the  parameter  X.  If  sample 
size  has  been  tabulated  as  function  of  lot  size  for  one  value  of  X;  X  *  1  say, 
we  may  therefore  use  the  same  tabic  to  find  the  sample  size  corresponding  to  lot 
size  N  and  any  X  by  using  N*  "  NX  as  argument. 


j,  Bayesian  single  sampling  plans. 

The  Bayesian  solution  consists  of  determining  the  value  of  (n,c)  minimizing 
R(N,n,c)  and  using  this  sampling  plan  if  min  R  is  less  than  the  costs  of 
accepting  or  rejecting  all  lots  without  inspection,  A  necessary  condition  for  a 
sampling  plan  to  exist  is  that  >  0  and  72  >  0,  i,c.  Pi  <  Pr  <  P2* 

Values  of  (n,c)  minimizing  R  must  satisfy  the  two  inequalities 
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^^R(N,n,c-l)  i  0  <  A^R(N;n,c),  0  S  c  S  n, 

and 

A^R(N,n-l,c)  S  0  <  A^R(N,n,c),  c  S  n  S  N, 

A  denoting  the  usual  forward  difference  operator.  Noting  that  A^B(c,n,p)  » 
b(c  +  l,n,p)  and  A^B(c,n,p)  -  -pb(o,n,p)  we  find  from  (19) 

A^R(N,n,c)  -  (Nn)(-7jb(c+l,n,pp+  72*’(c+l,n,P2)) 

and 

A^R(N,n,c)  -  1  -(7^Q(Pj)-*-72P(P2)^  (N-n-l)(^^pjb(c,n,p^)-y2P2‘’(^^">P2^^ 

Solving  the  inequalities  with  respect  to  n  and  N  we  find  that  a  Bayesian  sampling 
plan  must  satisfy  the  two  inequalities 


a  +  pc  ?  n  <  c:  +  p(c+l) 

(36) 

and 

F(n-l,c)  «  N  <  F(n,c) 

(37) 

where 

72 

a  •  (In— )/(ln“), 

'1  ^*2 

P2‘’l  ^1 

P  »  (In  1  (In  -f)  , 

qzPi  q2 

(30) 

(39) 

l-7,+rjB(c,n,pp-  72B(c,n,P2) 

F(n,c)  -  n  +  1  +  ; - r- - TT - . 

“7iPjb(c,n,Pi)  +  72P2‘^(^^"^P2^ 

For  two  plans  (n^^c^)  and  (02,02),  Cj  <  C2  say,  satisfying  (36)  and  having  over¬ 
lapping  N- intervals  according  to  (37)  the  cost  functions  must  be  compared.  Solving 
the  equation  R(N,nj,Cj)  -  R(N, 02,02)  for  N  we  get 

(n2-nj)(l-7p-!-  027(02,02)-  n^7(nj,Cj) 

^12  ■  7(02,02)  -  7(nj,Cj) 

where 

7(n,c)  =  7jB(c,n,Pj)  -  72^(0, n,P2). 

Since  u  for  given  (n,c)  is  an  increasing  linear  function  of  N  we  have  that 
R(N,nj^,Cj^)  >  R(N,n2,C2)  according  as  N  >  N^2» 

It  will  be  noted  that  the  simplicity  of  the  solution  depends  essentially  on 
R(N,n,c)  being  a  linear  function  of  N  for  given  (n,c). 

The  above  solution  has  previously  been  given  in  |^8  ]  and  a  rather  complete  tabulation 
of  Bayesian  sampling  plans  has  been  provided  in  |ll]. 
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The  exact  solution  given  above  does  not  disclose  the  structure  of  the  relation¬ 
ships  between  N^n,  and  c.  We  shall  therefore  derive  an  asymptotic  solution. 

Setting  c/n  ■  Pq  +  0  for  n  — >  we  shell  first  find  and  e  and 

afterwards  determine  the  relation  between  N  and  n  by  means  of  theorem  4. 


Using  (24)  and  the  expansion  of  b(c^n^p)  we  get  for  Pq  P 


B(c,n,p) 


‘*oP 


-1, 


p-PqI  .^Po^^o 


exp  (-ntp(pQ4€,p)  +  0(c)  +  p(n  )) 


and  the  same  expression  for  1  -  B(c,n,p)  for  p^  >  p.  For  P^  ^  Pq  ^  P2  therefore 
have  the  following  asymptotic  expansion  for  the  cost  function 


where 


and 


-1/2  ^  -1 
R  •  n  +(N-n)n  1  X  cxp(-n(p(p.,p.  )-ne(pJ  +  0(c)  +  0(n  )) 

i«l  ^ 


^i  “  Vi^  >  IPo'Pi' ^"Po^^o 

-  <p'(Pq,P^)  " 


(40) 

(41) 

(42) 


By  means  of  this  expression  and  theorem  4  we  shall  prove 
Theorem  5.  For  the  Bayesian  single  sampling  plans  we  have 


where 


.2 

nPo  +  +  a^yn  +  0(n  ) 

(43) 

"•l  P2‘’l 

(In  -i)  /  (In  ~)  , 

(44) 

^2  Pl^2 

‘  1  *  *  1  1  ^  "2  .  ^2  , 
‘2"  w  ^4“  ^ 


Po^i 


Po'i2 


(Pq-Pi)^  (P2-lV^  " 


(45) 


Further 


In(N-n) 


min  R  -  n  +  ~(1  -  — )  +  0(n  ^), 

^0  - 

'P^n+ ^In  n+ n^tpj  -  ln((j)^A.jO Vl-ipJ))  +('^^j'''{+3i’ +0(n”^) 


7  (1  -  ^  +  0(n  ")) 

'  2(p^n  —  N-n 


P(P2) 


1  1-2  1 

7  (I  -  T-^  +  0(n  ))  77^ 

7  '  ’  n  —  '  N-n 


V2^'  2cp^ 


(43) 

(4?) 

(40) 

(4S) 


0 


and 


where 


.  '’0^  . 

■^0  ■  Po  "0  ‘  ■ 

1  or  2  , 

(50) 

and 

2. 

“2+“2  ‘■Po''o  “2  , 

Po**i 

‘  *  ^Poio  ‘^Po’o  ’  Pi-Po 

(Pi-Po)' 

(51) 

Proof. 

For  Pj  Pq  P2  follows  from  (40)  that  R/N 

— >  0  for  N  — >  08, 

n  ”■>  00, 

nnd  n/N  — 0.  For  <  Pj^  we  have  that  Q(Pj)  — >  1  and  P(P2)  0  ao  that 

R/N  — >  0.  For  p^  -  p^  we  may  uae  the  normal  approximation  which  gives  that 
n(pp  — >a  ,  3  <  u  <  1,  so  that  R/N  —— >  >  0.  The  optimum  value  of  p^  satisfies 

consequently  the  inequality  Pj  Pq  P2* 

It  is  easy  to  sec  (indirectly)  that  the  optimum  value  of  p^  must  satisfy  the 
equation  ^(p^^Pj^)  =  fp(PQjP2)  which  leads  to  (44)  because  any  ocher  choice  of  p^ 
will  make  one  of  the  exponential  terms  in  (40)  larger  than  cxpC-ntp^^). 

To  find  a  first  approximation  to  c  we  put  dR/dc  =  0  which  gives 

2 

r  X.cp^  c  ^  0. 

i=l  ^  ^ 

Solving  for  a^  ne  we  get  the  first  part  of  (45),  i.e.  c  np^  f  82  +  o(l)  as 
could  be  expected  from  (3;j)  since  p^  *  1/P. 

To  geL  one  more  term  we  write  c/n  ^  p^  +  (a24c)/n,  e  — 0,  and  repeat  the 
procedure  ^o  find  the  new  c.  Using  the  same  method  as  in  section  3  we  may  develop 
(40)  into 

2 

-1/2  ^  -2 
R  n  ■l-(N  n)n  cxp{  •  nc,  -a-(pf  -  (nr  )/n  -{•  0(c/n)4O(n  )) 

i.-l  ^  2  i  I  1 

where  has  been  defined  in  (51).  From  dR/dr  »  0  we  find  the  equation 
ncff|+  nC(p'+P2  •  Solving  for  a^^  n£  we  get  the  second  part  of  (45), 

The  remaining  part  of  the  theorem  is  a  direct  consequence  of  theorem  4  for 

"^0  '■  ''o'  ^2  *  ^2’1  ■  ^ 

\  '\'^2  ^2'  ^  ^/(-vjtpp. 

It  will  be  noted  that  asymptotically  the  ratio  between  the  consumer's  and  the 
producer's  risk  is  constant  for  the  Bayesian  plans,  viz. 

^(P2)/Q(Pi) 


cf.  Theorem 


0 


From  (36)  one  could  hope  that  n  *  a  +  0(c  +■  and  consequently  that  *  "Pq^’ 
In  [ll]  it  has  been  demonstrated  numerically  that  this  approximation  is  rather 
good.  We  shall  now  derive  the  limiting  values  of  a^  and  for  p^  — 0  and 
^  ~  result  to  discuss  the  approximation  above. 

Proceeding  as  in  section  A  we  find 

ln(7^/72)  I 

^2  ^  In  r  '  ^  In  r  g2l'^((82''‘^^^82^ 


and 

Po<^4  — -(8j  +  82X^2  ■*"  h  ‘ 

The  last  term  in  lim  a^  decreases  from  0.50  to  0.A7  as  r  increases  from  1  to  2G, 
i.e.  the  last  term  is  practically  equal  to  1/2. 


For  small  p^^  and  we  therefore  have  approximately 


nPo  + 


Iniy^/y^) 

Tn~ 


2  np^ln  r  ''' 


The  Traction  of  the  average  decision  loss  which  is "due  to"  chc  consumer's  risk 

equals  (p|/6’.  This  fraction  Lends  to  (In(r-l)  In  In  r)/ln  r  Tor  — >  0 

(r  -  P2^^1  increases  from  0.50  to  0.62  as  r  increases  from  1  to  20. 


7.  Rescricted  Bayesian  sampling  plans. 

One  of  the  objections  against  the  B'yesian  solution  is  that  it  does  not  always 
lead  to  a  sampling  plan,  particularly  for  small  lots.  In  such  cases  a  running  check 
on  the  assumptions  regarding  the  prior  distribution  is  lacking  and,  if  quality 
deteriorates,  the  delay  before  appropriate  measures  can  be  taken  may  be  excessive. 

There  may  also  bo  cases  where  a  producer,  say,  inspecting  his  own  goods  sets  an 
upper  limit  for  the  probability  of  passing  bad  lots. 

It  is  therefore  useful  to  study  restricted  Bayesian  solutions  derived  by  minimizing 
average  cos^^o  under  a  suitably  chosen  restriction.  Such  restrictions  may  be  of  an 
economical,  technical,  or  statistical  nature.  We  shall  here,  however,  cnly  consider 
restrictions  on  the  operating  characteristic,  i.e.  restrictions  which  arc  independent 
of  the  weights  in  the  prior  distribution  and  the  cost  functions. 

The  first  restriction  of  that  kind  was  introduced  by  Dodge  and  Romig  [  A  J  in  their 
LTPD  system  of  sampling  plans.  For  the  present  model  it  consists  in  specifying 
P(P2)  3;  where  6  customarily  is  chosen  as  0.10.  Correspondingly,  one  may  choose 

to  specify  P(pj^)  =  a  with  'i  -  C.Oj,  say. 


19  - 


It  follows  from  the  pTopertieR  of  the  Bayesian  sanpllng  plans  that,  at  least 
for  large  lots,  it  will  be  uneconomical  to  specify  a  fixed  risk  for  the  consumer 
or  the  producer.  Instead,  the  risk  should  be  chosen  as  r.  decreasing  function  of 
lot  size,  for  example  as  P(p2)  ■  Pq  for  N  <  and  P(P2)  ■  ^  ^ 

similarly  as  Q(Pj)  •  for  N  5  and  Q(Pj)  ■  a^N^/H  for  N  >  N^,  which 
asymptotically  corresponds  to  (AO)  and  (A9), 

Another  possibility  is  to  specify  P(Pq)  *  1/2  for  Pj  Pq  ^  P2  '^I'lch  will  lead 
to  decreasing  risks  both  for  the  producer  and  Che  consumer.  If  p^  is  chosen 
as  in  (AA)  both  risks  will  be  0(1/N),  otherwise  they  will  tend  to  zero  with 
different  rates  of  convergence. 

Finally  one  may  specify  the  ratio  between  the  two  risks,  i.e.  Q(Pj)  *  pP(P2)> 

which  also  results  in  decreasing  risks  with  increasing  lot  size.  The  corresponding 
Bayesian  result  is  given  by  (52). 

It  will  be  seen  that  it  is  possible  (inspired  by  the  asymptotic  properties  of  the 
Bayesian  sampling  plans)  to  introduce  restrictions  on  the  operating  character^ at jr 
which  will  change  the  Bayesian  solution  for  small  lots  in  a  dasired  direction  and  a 
the  same  time  preserve  the  properties  of  the  Bayesian  solution  for  large  lots. 

It  should  be  noted  chat  restricted  Bayesian  plans  exist  also  if  one  of  the 
parameters  7^  and  ^2  1^  zero  which  means  Chat  the  double  binomial  prior 
distribution  degenerates  Co  a  single  binomial.  The  other  quality  level  is  then 
introduced  through  the  restriction  with  the  purpose  to  avoid  too  heavy  losses  if 
the  prior  distribution  changes.  In  sections  8  -  11  we  shall  discuss  the  properties 
of  the  various  types  of  restricted  Bayesian  sampling  plans. 


8.  .Minimum  average  costs  for  fixed  consumer's  or  producer's  risk. 

We  shall  limit  the  mathematical  discussion  to  Che  cesc  with  a  fixed  consumer's 
risk  since  the  two  cases  are  completely  analogous. 

Inserting  P(P2)  =  3  Into  (19)  we  get 

R(N,n,c)  '  n(l  -  '  n)7j^Q(p^)  +  NyjP 

-  (1  -  72P)Ri  +  NyjP  (53) 

where 

Rj  -  n  +  (N  n)7Q(p. ),  7  '  7^/(1  -yjP).  (5A) 

The  problem  is  to  minimize  R  or  equivalently  under  Che  restriction 


*'1 
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P(p2)  ^  BCc^n^p^)  P. 

Since  (55)  defines  a  relation  between  n  and  c,  n  n^  say,  (which  is  unique 
if  n  io  determined  as  the  smallest  integer  satisfying  B(C;n,P2)  5  P)  we  may 
write  Rj^  as  a  function  of  c 

R^(c)  -  n^  +  (N  -  n^)7(l  -  B^) 

where  B^  =  B(Cjn^,Pj^).  The  value  of  c  minimizing  is  determined  from  the 
inequality  /^^(c  1)  5  0  <  .  nj^(c).  As 

-^Rj(c)  -  (1  -  7  +  7B^)/n^  -  7(N  -  "c+l^^'^c 

we  define  the  auxiliary  function 


N  =  n  +  (1/7  -  1  +  B  )/'n  /^B 
c  c+1  '  '  c  c  c 


so  that  (n  ,c)  is  the  optimum  plan  for  lot  size  NifN  ,§N<N,/^B  ,>0 

c  r  r  J.  ^  c'  C  1 

and  /'^  >  C. 

c 

These  fonnulas  are  well  suited  for  a  systematic  tabulation  since  n^  and 
arc  easily  found  for  c  C,  I . 

The  above  solution  has  previously  been  discussed  and  tabulated  in  ^12^  •  A 
similar  approach  has  been  used  in  [s] ,  and  by  Dodge  and  Romig  .  It  will 
be  noted  that  the  Dodge  Romig  LTPD  plans  are  obtained  for  7j^  -  I  and  72  0  with 

the  modification  that  Dodge  and  Romig  use  the  hypergeometr ic  distribution 
Instead  of  the  binomial  in  defining  the  restriction. 

The  asymptotic  solution  is  given  in 


Theorem  G.  The  relation  between  acceptance  number  and  sample  size  is  given 

A 


.  ^  -i/2  -5/2. 

c/n  -  p„  -f-  I  a.n  +  0(n  ) 

/.  .it  ” 

1  1 


where 


,  ,1/2 

=  -u(p,q,)  , 


2 '*2 


a^- -~(l-6p2q2)(P2q2^  ^^^(u^-3u)  +3^6(q2'P2^  (P2q2^  ^  > 


(5:>) 


a 


‘^2’P2  I  I 


A  P2q2 


^(l-12p2q2)(u^-6«^+3)-  2) 


and  u  denotes  the  of  the  standardized  normal  distribution. 


The  relation  between  sample  size  and  lot  size  is  given  by  (33)  or  (35)^  and 
min 


2 

R.  +  (57) 

0  o'  4a„n  ^ 


and 


Q(Pl) 


iiii?  r. . 

’'1=0  V 


2-nVn  /  2 
0’ 


(5G) 


vhc 


re  «  (p(P2/Pi)^  \  for  i  j  1, 


and  X  =  (7q2Pi)/((P2*Pi)\.'2ffp2q2)- 


Proof .  The  asymptotic  solution  to  (55)  may  be  found  from  the  Fishcr-Corrlah  [ij 

2  2 

expansion  which  leads  to  (56)  since  k  /a  ■  q  •  P#  =  1  -  6pq,  and 

2  ^ 

(q-p) (l-12pq)  for  the  binomial.  The  u-polynomials  have  been  tabulated 
as  functions  of  p  in  [oj. 

The  remaining  part  of  the  proof  consists  of  a  direct  application  of  theorems  1 
and  4, 


Analogous  results  for  the  case  with  ■  1,  ^2  "  ^  and  Poisson  probabilities 
have  been  given  in  [9]  which  also  contains  an  evaluation  of  the  accuracy  of  thj 
asymptotic  solution.  In  the  present  paper  the  expansion  has,  however,  been 
carried  two  terms  further  than  in  [sj  and  [l2J, 

Since  the  producer's  risk  tends  (exponentially)  to  zero  and  the  consumer's  rick 
is  fixed  the  average  decision  loss  tends  to  a  constant  plus  72^^*  This  last 
term  will  for  large  N  become  dominating  also  as  compared  to  the  costs  of 
sampling  (which  are  of  order  In  N)  and  for  this  reason  plans  with  fixed 
consumer's  or  producer's  risk  arc  uneconomical  for  large  N, 

An  Important  result  is  that  asymptotically  the  sampling  plan  depends  only  on 
the  product  of  lot  size  and  cost  parameter  so  that  the  plan  for  lot  size  N  and 
cost  parameter  7  equals  the  plan  for  lot  size  Ny  and  cost  parameter  1.  Since 
this  property  holds  with  good  approximation  also  for  small  N,  tabulation  of 
such  plans  may  be  limited  to  7  ^  1  and  7=5,  say,  for  LTPD  plans,  and  to 
7  -  I  and  7  '  0.2  for  AQL  plans. 

The  present  LTPD  model  generalizes  the  Dodge-Romig  system  of  LTPD  plans  in 
two  respects:  (1)  the  single  binomial  prior  distribution  has  been  replaced 
by  a  double  binomial,  (2)  the  sim./le  cost  function  I  -  n  +  (N-n)Q(p^),  which 
pertains  to  rectifying  inspection 


with  the  same  costs  for  sampling  Inspection  and  sorting^  I  as  been  replaced 
by  R  which  allows  a  much  broader  interpretation  including  rectifying  as  well 
as  non-rectifying  inspection. 

In  case  w^*  0,  i.e.  0,  the  quality  level  p^  is  introduced  through  the 
restriction  alone,  which  means  that  P2  has  to  be  deteminod  from  technical 
and  economical  considerations  and  not  from  the  prior  distribution. 

A  detailed  discussion  of  both  LTPD  and  AQL  plans  and  a  comparison  with  other 
systems  has  been  given  in  [12]. 


9.  Minimum  averap.c  costs  for  P(Pq)  ■  1/2. 

The  restriction 

P(Pq)  »  B(c,n,pQ)  -  1/2,  p^<  Pq<  P2, 

defines  a  relation  between  n  and  c,  n  =»  n^  say.  Proceeding  as  in  section  C  we 
find  that  the  plan  (n^,c)  is  optimum  for  1  ^  ^c 


N  - 

c 


n 


c*!-i 


(1-7^Q(Pj) 


72P(P2))^*'c/(^j^(Pi)  + 


72^(P2))  ^ 


Q(p^)»  1  -  IJ(c,n^,pp  and  P(p2)  =  B(c,n^,P2). 


In  principle  such  IQL  plans  may  be  defined  for  any  value  of  p^.  If  no  (technical) 
reasons  exist  for  choosing  a  .specific  va’uc  ol  it  seems  reasonable  to  determine 
Pq  such  that  R  Is  minim ’.zed .  From  (26)  i  Collows  that  one  of  the  risks  will 
be  Infinitely  snuall  as  conpared  to  the  ot.her  for  n  —>  t»  depending  on  which  of 
the  two  coeffictcnco  cp(pQ,p^)  and  (f(Pq,P2)  is  the  larger.  If  (p(pQ,P2)  ^  <^'(Pq>Pj) 
we  find  that  min  R  ~  r  l/.^(pQ,pp, 

-9(Pr,,P7)/(f(Pn,P,) 

Q(pp  -  0(1/M),  and  P(p2)  «  0(N  ^  ^  ). 

Since  n  ~  (In  N)/(p(pQ,pp  it  follows  chat  min  R  will  be  minimized  with  rcsiiect 
to  Pq  for  g(pQ,pp  being  as  large  as  possible,  l.c.  for  (p(Pq,p^)  =  ‘p(Pq>P2^' 
which  leads  to  as  defined  by  (^A) .  We  shall  therefore  only  study  IQL  plans 
for  this  value  of  p^. 

Theorem  7 .  The  relation  between  acceptance  number  and  sample  size  is  given  by 


2_2o  ,  ,,  -2. 

'"’o'  3  ■  32''OnpQq|j  *  2(n  )  • 


(59) 


The  relation  between  sample  size  and  lot  size  Is  given  by 


w 
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and 


n,ln  R  -  n  +  ^(1  -  j;|^)  +  0(n‘*), 

h  -“z”!  / ,  .  1 


1  +  i<b,  -  b, ,  -  ri-)  +  0(n‘*)^  ^ 

n  4  41  2(fi'  “  '  N-n 


Q(p,)  ■  — c  1  I 

(Pq\7,  V 


where  a^-  -(2-Pq)/3,  ^  ^2*^  * 


r 


•a  (p!  -a 


(61) 


(62) 


Proof ■  The  first  patL  of  uhc  theorem  follows  from  (56)  for  3  -  1/2,  i.c.  u  ■  0. 

2 

Let  c/n  “  Pq+  a2/n  +  a^/n  where  and  are  defined  by  (59).  Then  from  theorem  1 
we  have 

-1/2  -2 
R  n  +  (N-n)n  '^Xe  (I  +  0(n  )) 


with  the  definitions  of  X  and  b^  as  given  above.  The  remaining  part  of  the 
proof  follows  from  Lheorcu  4. 

Since  a^  is  small  we  have  practically  c  a  npQ-2/3  and  consequently  PCp^)^  Q(pj) 
according  to  theorem  2. 

The  IQL  plans  are  the  only  plans  with  a  fixed  risk  which  lead  to  a  relationship 
between  c  and  n  of  the  same  type  as  the  one  for  the  Bayesian  plans.  The 

coefficients  82  and  a^  arc,  however,  different.  As  a  result  we  find  in  the 

IQL  system  P(P2)/Q(Pj)  ^  1  whereas  this  rado  in  the  Bayesian  system  equals 

Because  of  the  relation  Q(Pj^)  *  P(p2)  have 

R  a  n  +  (N-n)(7j+  72)Q(Pp 

so  that  we  may  use  chc  IQL  plans  defined  by  (Pq>  y-^  approximation 

to  the  IQL  plans  defined  by  (Pj,P2>7j>72^  which  makes  a  tabulation  much  easier, 

sec  [12].  Plans  for  7.+  7-^  1  may  be  found  with  good  approximation  from  a  table 

for  7^+  72*  1  by  using  N  “  N(7j+  as  argument. 

The  IQL  plans  considered  here  are  generalizations  of  the  plans  discussed  by 
VJcibull  [17]  and  tabulsLcd  by  Markback  [iS]  in  the  same  manner  as  the  LTPD 
plans  are  generalizations  of  the  Dodge-Romig  plans. 


10.  Minimum  average  costs  for  decreasing  consumer's  or  producer's  risk. 

The  restriction  imposed  is  of  the  form  P(P2)  ■  P(N),  say,  where  0(N)  is  a 
given  decreasing  function  of  N  tending  to  zero. 


•»  f 
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We  shall  only  treat  the  case  P(P2)  ■  Pq  for  N  S  and  PCp^)  " 

for  N  >  Nq.  For  a  specified  decreasing  producer's  risk  we  get  analogous 

results . 


For  N  S  Nq  the  theory  has  already  been  given  in  section  8.  For  N  >  the 
procedure  for  obtaining  the  exact  solution  is  analogous  to  the  one  derived 
in  section  8  with  the  (numerical)  complication  that  n,a3  defined  Ly  the 
restriction,  becomes  a  function  of  both  c  and  N,  Therefore  the  auxiliary 
function  N  has  to  be  determined  by  iteration  for  each  c. 

The  asymptotic  solution  is  given  in 

Theorem  8.  The  rela.ion  between  acceptance  number  and  sample  size  is 

-2. 


c  •  npQ+  a^d-  3^/n  +  0(n  )  where 


1 


1 


1 

a^^  ~  In  -  and  a  =  —  (P  -  3,+  7^). 

2  6’  XjC-gp  '•  B'  >■  ‘  ^’0 


(64) 


Po>  Vq,  and  P^  being  defined  in  (44),  (50),  and  (51),  reopecuivcly . 
The  relation  between  sample  size  and  lot  size  is  given  by 

1  ^^2  1-2 
In  N  ■  cp^n  +  -  In  n  ■¥■  In  P2^n  ^ 


and 


^“^2^  1  -2 
min  R  -  n  d-  Py^  +  d  *  +  0(n  )  . 


(65) 


(66) 


Proof.  Writing  c/n  -  Pq^-  c  it  follows  as  in  section  6  that  Pj^  <  Po<  P2 


From 


B(c 


\  .1/2  -n9(P(,,Pj)-nc<t’  ^ 

,n,P2)  -—no  .  -Jp  , 


disregarding  terns  of  "higher  order",  wc  find 


ci(pO»P2^"  +  jln  n  +  neep'  ln(372/^2^  * 


(67) 


For  given  values  of  p^  and  c  this  gives  a  relation  between  n  and  N.  Il  remains 
to  determine  p^  and  s  so  that  R  is  minimized. 


Writing 


R  -  n  +  (N-n)72P(P2)  ((7jQ(Pi))/(72P(P2)) 


/\  -nb  -nr5' 

R  =  n  -I-  (I  -  ^)P7  ’  —  c  +1 

N  2  ^^2 


wc  find 
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From  (j7)  wc  have  n  (In  N)/(p(p  ,p-)  eo  that  p  enters  into  R  through  n  and 

0  c 

expC-nb^).  To  determine  the  optimum  value  of  p^  consider  first  6^  0  which 

gives  p  as  determined  by  ('»A).  As  shown  below  we  then  determine  c  so  that 

”  n6.  H 

10+  0(1).  For  8q  <  0  wc  get  e  ^  -  0(ir),  u  >  0,  which  gives  a  larger  R 

thru  for  6^  0.  For  8^  >  0  wc  find  R  n  +0(1)  but  since  tp(PQ,P2)  <  Vq  the 

corresponding  R  will  be  larger  chan  for  6^  -  0.  Consequently  we  have  6^  C. 

To  determine  c  wc  solve  the  equation  dR/dr  -  0  taking  into  account  that  n 
according  to  (j7)  depends  on  c.  Using  dn/^c  ^  solving  for  a^  nr 

wc  get  the  result  given  in  (JiA). 

To  get  one  more  term  we  introduce  c/n  *  Pq  (^2  ^  repeat  the  procedure 

to  determine  the  new  c.  From 

(f^n  +  -|ln  n  +  ^  ^  I*'  N 

and 

n  /^i  a  8'-(j5  -p.+nrBO/n  \ 

wc  get  dn/or  =  (-  ~  and  solvin’  the  equation  dR/dc  ’  0  with 

respect  to  nc  we  fin'’  as  given  in  (£>A). 

The  remaining  part  of  the  theorem  follows  by  inserting  the  values  of  a^  and. 
in  R. 

This  system  of  plans  has  thos  the  same  asymptotic  properties  as  the  Bayesian  plans. 
The  essential  difference  between  the  two  systems  arises  from  the  different 
constant  terms  in  the  asymptotic  relations  between  c  and  n  and  between  n  and  N 
because  i.aecc  constants  depend,  cu  the  (arbitrary)  parameter  (3. 

V’c  may  minimize  min  R  with  respect  to  3.  Taking  into  account  that  n,  a^,  a^,  and 
^2  are  functions  of  3  we  find  the  optimum  value  is  3  '■  '<’1^  ^  might  be 

expected  from  the  Bayesian  solution  and  furthermore  min  R  *  n  +  ^)* 

11.  Minimum  average  costs  foi  a  fixed  ratio  of  the  consumer's  to  the  producer's  ri.sk . 

The  restriction  is  given  as  P(P2)  ~  pQ(Pj)  which  naturally  may  be  modified  so 
that  for  N  5  Nq  it  is  further  required  that  P(P2)  =  3>  and  that  P(P2)  S  3 

for  N  i  Nq. 

The  restriction  defines  n  os  a  function  of  c  and  the  exact  solution  may  therefore 
be  derived  by  the  same  method  as  used  in  section  C.  Since  P(Pj)  >  P(P2)  it  is 
necessary  for  a  solution  to  exist  that  P(P2)  <  p/(l+p). 


Tho  flsympCoClc  solucion  Is  given  by 


io  - 


Theorem  9,  The  relation  between  occcptance  number  ood  sample  sice  is 

■  2 

c  ^  np^  +  -f  a^/n  +  0(n  )  where 


1  , 
In 


5<  -  ' 


p^  and  a^  being  defined  in  (44)  nnd  (30),  respectively. 
The  relation  between  sample  size  and  lot  size  is  given  by 


I  xl 


-2. 


In(N-n)  -  -In  n+  -  ln(X(,iQ)+  (a^ipj  +0^  -  2ip~^n  ^ 


and 


min  R 


n  +  —  (I 


0 


(6C) 

(-9) 

(7D) 


where  X  *  defined  in  (51). 

Proof.  The  proof  is  similar  to  the  one  in  section  6  making  use  of  theorems  3  nnd 

4. 


Minimizing  min  R  with  respect  to  p  gives  0  =  (y^tpp/ (-729^  which  leads  back  to 
the  Bayesian  solution. 

For  p  1  we  find  for  small  and  that  a^  a  -2/3  which  means  that  we  get 
approximately  the  IQL  plans  discussed  in  section  9. 


12,  Sampling  plans  defined  by  two  risks. 

Suppose  now  that  the  weights  of  the  prior  distribution  and  the  costs  are  unknown. 

Uu  then  ask  the  question:  Is  it  possible  from  knowledge  of  (Pj^,P2)  alone  to 
define  systems  of  sampling  plans  having  similar  properties  os  the  Bayesian  and 
the  restricted  Bayesian  plans? 

By  means  of  the  asymptotic  properties  derived  in  the  previous  sections  we  may 
construct  three  systems  with  the  required  properties,  each  of  thonbeing  defined 
by  specifying  two  risks: 

A.  Decreasing,  producer's  and  consumer's  risk,  i.e.  Q(pp  *  a/N  and  P(P2)  “  P/N, 
which  corresponds  to  the  Bayesian  plans  in  section  5  and  the  restricted 
Bayesian  plans  in  section  10  and  11. 

B.  LTPD  plans  with  decreasing  producer's  risk  (or  AQL  plans  with  decreasing 
consumer's  risk),  i.e.  Q(pp  =  c/N  nnd  P(P2)  ■  P  (or  Q(Pj)  ■  a  and 
P(P2)  •  0/N);  which  corresponds  to  the  plans  discussed  in  section  0. 
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IQL  plans  with  decreasing  producer*a  (or  contumcr^s)  risk,  l.c. 

PCPq)  ■  1/2  and  ©(p^)  ■  a/N  (or  P(P2)  “  3/N),  which  corresponds  to  the 
IQL  plans  discussed  in  section  9. 


It  may  be  necessary  to  modify  conditions  of  the  form  ©(p^  ■  a/N,  say  ,  to 
Q(Pl)  Oq  for  N  i  Np  and  0(pp  »  ^  ^  ^0' 

It  will  be  noted  that  the  requirement  min  R  has  been  replaced  by  conditions  as 
Q(Pj)  a/N  and/or  P(P2)  P/N* 


Wc  shall  give  a  short  discussion  of  these  three  systems.  Each  of  the  first  two 
contains  two  arbitrary  parameters,  a  and  p,  whereas  the  last  contains  only  one. 


A,  From  the  two  conditions  wc  find  that  B(c,n,P2)/(l  '  B(c,n,Pj^))  •»  (3/a,  which 

defines  a  relationship  between  n  and  c,  n  *  n^  say.  Defining  the  auxiliary  function 

N  •“  3/B(c,n  ,p-)  It  follows  that  (c,n  )  satisfies  the  conditions  for  N  ,  <  N  <  N 
cc2  c  c-1  c 

(If  we  Interpret  the  equalities  as  S).  It  is  therefore  rather  easy  to  tabulate 

thcffc  plans  when  n^  has  been  found.  Since  P(pj)  >  P(P2)  ft  is  necessary  for 

a  solution  to  exist  that  N  >  a  +  p. 


2  -3 

Asymptotically  wc  find  that  c/n  “  Pq  ^  where  p^,  a^,  and 

are  determined  by  means  of  theorem  3,  l.e. 


a 


2 


JL  , 

6^  aX2ri 


(71) 


The  relation  between  lot  size  and  sample  size  is  given  by  (55). 

The  corresponding  costs  arc 

R  ^  n  +  (7ja  +  72^)  ^  ^)  (72) 

which  is  minimized  with  respect  to  a  and  3  for  a  -  (-(pp/ ,>j(}>q6'  and 
compare  theorem  3. 


For  any  a  and  3  we  may  thus  construct  a  system  of  sampling  plans  with  the  same 
asymptotic  properties  os  the  Bayesian  plans.  Tabulation  of  the 'bptimun/' values 
of  a  and  3  for  typical  values  of  (Pj,p2>7j#72^  give  some  guidance  for  the 
choice  of  a  and  3. 


B.  The  condition  T^Oo)  "  P  defines  n  *•  n  which  together  with  the  other  condition 
—  '■  I  c 

gives 

N^  =  a/(l  -  B(c,n^,Pj^) )  .It  is  necessary  that  N>  a/(l  -  3)» 


Asymptotically  c/n  is  given  by  (55),  and  by  means  of  theorem  1  we  get 
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In  N  -  cp(p2/Pj)n  +  ^Inn-fl^+lnCaCp^-p^)  +  ...  (73) 

where  -  b^Cp^^p^). 

The  corresponding  costs  become 

R  -  (1  -  +  7^^^  +  0(n"^)  (7A) 

%rhich  IS  minimized  with  respect  to  a  for  a  =  (I  -  compare 

theorem  6. 


The  system  may  be  characterized  as  an  LTPD  system  with  the  consumer's  risl;  inversely 
proportional  to  lot  size.  It  is  much  easier  to  tabulate  than  the  correspondinj 
restricted  Bayesian  system. 


C,  The  exact  solution  is  analogous  to  the  one  under  B. 

2  -  3 

Asymptotically  we  have  c/n  “  Pq  +  ^2^*^  "*■  ‘’*4^^  "*■  where  the  cocf ricienir. 

a^  and  a^  are  given  by  (59).  From  theorem  I  we  find 


In  N 


1  -1-2 
(p^n  +  yin  n  +  a^cpj  +  InCay^/X^  -f  (p^^p^n  +  0(n  ). 


i'j) 


The  coats  become 

a  b' 

R  •=  n  +  7^0  +  ^ — -  c  (I  +(a^6'  +  3^ 


which  is  minimized  with  respect  to  a  for  a  ■ 


-  Oj)n  ^ )  V  0(n  ^) 

■^2'^1 

XjC  /(p^X7j,  compare 


th^cren 


(lb) 


7. 


The  system  is  an  IQL  system  with  the  consumer's  risk  inversely  proportional  to  lot 
size.  According  to  theorem  2  the  producer's  and  the  consumer's  rick  are  nearly 
equal. 

The  exact  solution  is  very  easy  to  tabulate  because  may  be  found  from  the 
asymptotic  solution  with  sufficient  accuracy  for  c  ?  1. 


D.  Fixed  risks.  For  the  sake  of  completeness  we  mention  the  system  with  fixo'’. 
risks,  i.c.  Q(Pj)  ®  a  and  P(P2)  =  3^  which  leads  to  a  sampling  plan  indepencent  0 
N  .and  gives 

R  -  (1  -  7^a  -  y^^)x\  +  (y^x  +  72p)N*  (77) 

For  systems  A  and  C  the  costs  consist  of  sampling  inspection  costs,  n  =  0(ln  F); 
plus  an  asymptotically  constant  average  decision  loss. 

For  system  B  the  sampling  inspection  costs  arc  0(ln  N)  and  the  da:lsion  los'’ 
consists  o^  a  constant  part  (from  the  decreasing  ric!0  and  r.  part  proportional 
to  N  (from  the  constant  risk). 
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For  system  D  sampling  Inspection  costs  arc  constant  and  the  decision  loss  is 
proportional  to  N. 

E,  PercentaRC  inspection.  The  plans  arc  defined  by  n  =•  nN,  u  being  a  (arbitrary) 
constant,  and  c  «  p^n  (or  c  -  p^n  +  +  ...). 

From  theorem  1  we  get 

R  ^  uN  +  0(Ne"^)  (78) 

so  tli.nt  we  have  the  interesting  result  that  for  the  given  model  plans  with 
fixed  risks  ajid  percentage  inspection  asymptotically  have  the  some 

costs  for  4  - 

In  the  first  ease  the  costs  arc  essentially  decision  losses  (because  the  cample 
is  too  small)  and  in  the  second  case  the  costs  are  essentially  inspection  costs 
(because  the  sample  is  too  large). 

For  1 ,  y2  ■  0.  5,  a  ”  0.05, and  p  =  0.10  we  get  u  0.10  which  is  of  the  order 

of  magnitude  as  previously  fouod  '‘n  practice. 

13.  Efficiency  and  robustness. 

In  a  previous  paper  {lo)  it  has  been  proposed  to  define  the  efficiency  of  a 
sampling  plan  as 

c(N,n,c)  =  RQ(N)/R(N,n,c)  (79) 

where  Rq(N)  denotes  the  costs  of  the  optimum  plan  and  R(N,i»,c)  denotes  the 
costs  of  the  plan  in  question. 

For  a  lot  containing  X  defectives  acceptance  without  inspection  is  cheaper  than 
rejection  W' thout  inspection  if  X  5  [Np^].  Classifying  all  lots  in  this  way  the 
average  costs  become 

[np.]  b 

K  (N)  =  T.  (NA»,  +  XA-)f..(X)  +  I  (NR|  +  XR,)f„(X). 

X-0  ^  ^  ^  X-[Np^]+l  ^  ^  ” 

By  means  of  theorem  1 

K  (N)/N  -  k  +  0(N"^^^  e'*^) 
m  m  “ 

and  K  (N)/N  — >  k  from  below, 
m  m 

It  would  be  more  correct  to  define  efficiency  as  the  ratio  of  costs  in  excess  of 
K  (N)  instead  of  Nk  .  The  difference  in  the  two  definitions  is,  however,  of 
importance  only  for  very  small  N.  The  definition  (79)  tends  to  underestimate 
the  efficiency. 
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Theorem  10.  Lee  a  system  of  sampling  plans  be  defined  by  ^hc  cwo  rc*a  '.our. 
c  -  npQ+  a2+  o(l)  and 

In(N-n)  ■  (p^n  y  In  n  +  k  +  o(l)  (3'. 


where  p 


0  V 


'  In  — 


In' 


and  <p 


‘O  ^0 

•  p^lu  —  +  q^ln  —  ,1-1 
0  ^0  0  q.  ' 


or  2.  Ti.cn 


c(N,n,c) 


— ^  (  "  X.  +  K--  K  -  1)  +o(n  S 
(PqH  C  0 


('A) 


K“a2<,^j^  K 1 2 

where  Kq  denotes  the  constant  term  of  (A7)  and 

Proof .  From  theorem  !3  wc  have  that  Rq(N)  '•  0^+  "o  deftne’  by 

(GO)  for  K  ®  Kq-  By  means  of  theorem  1  and  (GO)  we  find  that  R(N,n,  ')  •  p  X 
so  that 


1  \  '  X.  -1 

c  =  1  +  +  —  1  +  -  +  o(n  )  . 

n  n 

Comparison  of  (GO)  and  the  correspondins  equation  defining  n  lenur,  co 

- 1  ^ 
n^/n  ■  1  +  (k-k^)  /  +  o(n  )  which  corai<ietes  the  proof. 

The  theorem  shows  chat  the  restricted  Bayesian  plan.s,  apart  from  .he  one  1: 'Vi.., 
a  fixed  risk,  and  dio  corresponding  plans  based  on  two  risks  all  have  asyig.-ot  n 
efficiency  1.  The  reasons  for  this  impotwan-  result  arc  that  those  plans  only 
differ  with  respect  lo  the  constanL  terns  and  k. 

I.,  should  also  be  no-cd  ..hat  and  denend  on  (p^p^)  nnly,  udicrcis  a^  nr  ’  r 
depend  on  the  other  parameters  also.  This  means  chat  wrong  values  of 
and  the  cost  parameters  have  a  secondary  influence  on  the  c  f '"ic  ienty  t’h’ch 
uends  to  1  if  only  (p^,P2)  correct. 

Since  plans  having  the  consurier’s  and/or  the  producer’s  risk  fi'tcl,  i  J 
to  costs  of  order  M,  -he  efficiency  of  such  plans  tends  to  zero  ac(i:'.  '  /  N. 


Suppose  now  that  plans  have  been  constructed  from  wrong  values  o'  he  uarp:;''tcrs , 
(pJ.pJj  where  pj  =f  and  p*  [  P2 .  We  then  get 

Theorem  1 1 .  Let  a  system  of  sampling  plans  be  defined  by  the  two  rcla-iens 

*  « 

c  ■  npQ  +  a2  +  0(1)  and 

In(N-n)  ■  cp^n  +  jlnn+K+o(l).  <’32) 

Then  *  * 

^*0  '  1  •  1  ^"*0  '  - 1 
c(N,n,c)  -  —  1  -  k.-k:  -1  +  -  In  —  +  o(n  )  Ul) 

'•"o  (p'n  °  ^ 

*  H  -f- 

if  and  only  if  Pl<  P2<  P2 .  Otherwise  c  ”•>  0  as  N  In  N  for  0  <  c  *•  1 . 


I 
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Proof .  If  Pj<  Pq  <  p^  if  follows  from  theorem  1  that 

R  n  +  (N-n)n 

n  +  rx^fc 

#  ##  #4#  ##  # 

where  ‘P(Pq,P^)-^'’(Pq,P^)  since  Vq=  tp(pQ;Pj)  =  <P(p^j,P2)  •  Pq  ^  (P^;P2) 
.hen  R  contains  a  tenn  proportional  to  (H-n)  beside  terns  as  above. 


Since 


1  Pi 

6,>  0  for  and  6.  >  0  for  ; 

^  1  ^  1 

IPi>  P2 


'  P2^  ^2 


P2<  Pi 


we  get  6^>  0  and  0  ior  (1)  p^<  p*<  p^<  (2)  p^<  p^<  p^<P2  . 

(3)  Pj<  P2<  pj^  P2;  &!>  0  snd  5^<  0  for  (/ )  Pi<Pj<  P2‘^  P2^  ^1*^  °  ° 

for  (5)  p^<  pj<  ?*<  p^;  0  and  62<  0  for  (0)  p*<  p^<  p^<  p^. 


In  eases  (2)  and  (3)  we  have  p^  ^  (Pi>p2^  wherefore  R  =  0(N)  ,  In  eases 
) 

h  <  0  gives 


we  nay  have  Pq  '  (Pi>l2)^  since  ac  leas^:  one  of  the  is  negaLive,  and 


exp(-nf))  -  cxp(-~  In  N  +  ...)  =  0(N  ),  0  <  c  <  1, 


0 


we  get  R  =  O(N^),  0  <  f  *^1,  including  ^hc  cases  with  p^  c^:  (Pi^P^)'  follows 

-hat  c  “>  0  as  (In  N)/!!*"  . 


In 


-n 


case  (1)  we  have  R  “  n  +  0(c  )  so  thac  c  '  (*^q+  1/c^/n.  Fron  (02)  and  the 


corresponding  expression  for  n^  we  find 


n 


"0  /  *  1  *^0  _1_ 

=  —  +  I  T  In  ~  +  0(n  ) 

\  fi  i  /•*.  rr> 


'0  2 


^0  ^  '"o" 


which  leads  to  (83) . 


1, 


It  is  also  clear  hat  the  eases  with  at  least  one  risk  fixed  give  c  "  0(M  In  N) 


The  general  conclusion  is  that  the  Bayesian  plans  (and  also  the  restricted 

«  « 

Bayesian  plans  with  decreasing  risks)  arc  rather  robust  if  only  Pl<  ?!<  P2* 
A  discussion  with  nuncri''al  examples  has  been  given  in  [11]. 


If  the  prior  distribution  of  p  is  continuous  we  get  asymptotically  that  n 
is  proportional  coyN.  This  ease  has  been  discussed  in  [ 10  ]  where  also  a 
comparison  with  IQL  plane  has  been  given. 


t*  »•»-.  * 
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14.  An  example. 

As  an  example  consiuer  a  case  wi^h  Pj"-*  0  Oi^  0.85,  and  0.05,  W2”  0.15. 
Let  the  costs  of  sampling  inspection  be  0.40  (economic  units)  per  item  in  the 
sample,  i.e.  Sj^«»  0.40  and  0,  the  cos:g  of  rejection  0.30  per  item  in 
the  remainder,  i.c.  R^-*  0.30  and  8^=  0,  and  chc  costs  of  ac:eptancc  per 
defective  item  10.00,  i.c.  A^-  0  and  A^^*  iO.OO.  It  follows  that  0.03, 

y ^  0.6296,  and  0.1111. 

In  the  table  we  have  compared  plans  from  C  systems  defined  as  follows: 


(1) .  Bayes.  Plano  minimizing  R(N,n,c). 

(2) .  IQL.  Min  R  for  P(Pq)  •=  1/2.  We  get  p^-  0.0250.  The  plans  have  been 

found  by  minimizing  R  ■  n  +  (N-n) (7^+072)  where 

,  ,  -2/3 

p  =  -  O.9C02,  and  7,+P7~=  0.74. 

S’'2\Pl‘»2/  '  2 

(3) .  LTPD.  Min  R  for  P(P2)  “  0.10.  7  =  0.64. 

(4) .  AQL.  Min  R  for  QCp^  =*  0.05.  7  ■»  0.11. 

(5) .  Fixed  risk.  Q(pj)  “  0.05  and  P(P2)  -  0.10. 

(6) .  Percentage  inspection,  p  »  0.057^+  0.1072=  0.04259, 

(7) .  Dodge.  The  AQL  system  with  57.  consumer's  risk  proposed  by  Dodge  in  [3] 

with  AQL  =  Pj^. 

(8) .  Mil-Std.  Military  Standard  105D  [16]  with  AQL  =  p^. 

The  Dodge-Romlg  systet.i  has  not  been  Included  in  the  comparisons  because  it 
gives  nearly  the  same  result  as  the  LTPD  system  with  7  •  0.64. 


For  each  of  7  lot  size.s  the  Bayesian  plan  and  the  corresponding  costs  have 
been  found  and  for  the  other  systems  the  efficiency  has  been  computed  from 
(79). 


Since  7^  is  considerabl/  larger  chan  72  will  be  expected  that  plans  w.th 
small  values  of  Q(pj^)  as  compared  to  P('r2)  to  be  preferred.  Systems  with 

a  fixed  producer's  ri.sk,  such  as  (4),  (5),  and  (7),  may  therefore  be  expected 
to  give  low  efficiencies  for  large  lots,  as  seen  in  the  table. 

It  is  Interestin,  co  compare  the  two  risks  for  Che  various  systems. 

The  table  reveals  what  price  must  be  paid  .:o  obtain  a  specified  degree  of 
protection.  It  may  in  some  circumstances  seem  reasonable  to  pay  such  a  price 
for  small  lots,  but  certainly  not  for  large  lott  where  the  protection  obtained 
in  terms  cf  the  two  risks  is  very  good  both  for  the  Bayesian  and  the  IQL 
system. 
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15,  MigccIIaneom  remarkn, 

A.  Gcncralizationf  of  the  model.  As  Indicated  in  section  2  the  model  may  Lc 
Cencrr.lircd  by  introducinjj  a  polynomial  cost  function  "•a  in  (12)  and  /  or 
a  more  general  prior  discribution.  The  asymptotic  results  hold  for  any  prior 
distribution  having  probability  zero  for  values  of  p  in  the  open  interval 
(Pl,P2),  where  Pj  Pj.  ^  P2^  assigning  finite  probabilities,  and  w^,  ;"y, 
Wj  +  S  1,  to  the  end-points  of  this  interval,  i.c.  the  prior  distribution 
may  be  arbitrary  outside  the  closed  interval  (Pj^,P2)^  2*-''-’ 

Another  possible  generalization  will  be  to  make  the  prior  distribution  a  function 
of  N,  see  (13),  for  example  to  use  a  double  binomial  with  P2/P^ 


B.  The  AOQL  system.  If  p^  <  p^^  <  p^,  Pj^  denoting  the  ACQL,  it  may  be  shown  ''s 
in  [13]  that  R  =»  n  +  +  0(N^  where  cp^  =  ^ 

For  Pj^  ^  p^  we  have  cp^  which  mccans  that  the  ACX)L  system  givcc  a  satisfactory 

result  within  the  given  framework  only  if  p^^  ?  p^.  A  generalization  of  the 
A(DQL  system  has  been  discussed  in  [12] . 


C.  The  relation  to  hypothesis  testing..  The  producer's  and  consumer's  ricks 
correspond  to  the  probabilities  of  errors  of  the  first  and  second  kind,  respect  wo ly, 
for  testing  the  hypothesis  p  Pj  .against  the  alternative  p  =  p2»  Lehmann  [K:  | 
has  suggested  .o  use  the  relationship  p  *=  '  a  as  a  rule  cf  thumb  for  obt.aining 
a  reasonable  balance  between  the  two  error  probabilities  (.  -  Q(p,)  and  p  ~  P(p^)) 
instead  of  just  using  a  standard  signif ic.ancc  le'vel  for  ..  It  is  interesting  to 
notice  that  the  Bayesian  solution  .'‘symptot  ical  ly  has  this  property  ,  see  (52), 
so  that  Lclimann's  rule  is  supported  not  only  from  the  minimax  point  of  view, 
as  me*ntioncd  by  himself,  but  also  from  the  Bayesian. 


Ackpow lodgement. 

My  thanks  arc  due  to  Professor  D.  Duguc,  I'Institut  dc  Stetistique  do 
I'Universitc  de  Paris,  who  invited  mo  to  give  three  lectures  on  the  theory  of 
sampling  inspection  in  April  1935,  The  present  paper  is  a  revised  version  cf 
the  lectures. 

My  thanks  are  due  to  Mr.  N.  Keiding  and  Mr.  P.  Thyregod  for  checking  the  .asymptotic 
expansions.  Mr.  Thyregod  .also  carried  out  the  computations  for  the  example. 


35 


References. 

1.  Blackwell,  D.  anrf  Hodr^cs,  J.L.  (1959).  The  probability  in  the  extreme 
tail  of  a  convolution.  Ann.  Math.  Statist.  30, i 1 13- 1120. 

2.  Brockwcl  1.  I*. J.  (196^«).  An  asymptotic  expansion  for  the  tail  of  a  binomial 
distribution  and  its  appliction  in  queueing,  theory.  J.  />ppl.  Prob.  1,  161-.j7. 

3.  Dodr'o,  H.F.  (19(j3).  A  general  procedure  for  sampling]  inspection  by  attributes  - 
base!  on  the  /iQL  concept.  /^SQC  /.nn.  Conv.  Trans.  1963,  7-19. 

I ,  DouKC,  H.F.  and  Romir.,  H.G.  (1929).  /.  method  of  sampling  inspection.  Bell.  Syst. 
Tech.  J.  n,  ')13-31.  Reprinted  in  Dodge,  H.F.  and  Romig  H.G.  (1959).  Samplin;: 
inspection  tables.  Wiley,  New  York. 

5.  Dod.  e,  H.F.  and  Romip, ,  H.G.  (1941).  Single  ianpling  and  double  campling 
inspection  tables.  Bell.  Syst.  Tech.  J.  20,  l-l.  Reprinted  in  Dodge.  H.F.  and 
Romig,  H.G.  (1959)  Sampling  j.n9pcction  tables.  V'ilcy,  Mew  York. 

6.  Fisher.  K.  A.  and  Cornish,  E.A.  (19'j0).  The  percentile  points  of  distributions 
having  known  cumulants.  Tochnometr ics  2,  209-225. 

Guthrie,  D,  and  Johns  M.V.  (1959).  Bayes  acceptance  campling  procedures  for 
large  lots.  Ann.  Math.  Statist.  30,  896-925. 

C.  Hald,  A.  (1900).  Tl.e  compound  hypcrgeometric  distribution  and  a  system  of 
single  sampling  inspection  plans  based  on  prior  distributions  and  costs. 
Technometr ics  2,  275-390, 

9.  Hal  1 ,  A.  (1962).  Some  limit  theorems  for  the  Dodge-Romig  LTPD  single  sampling 
inspection  plans.  Tcchnonctrics  4,  497-513, 

10.  H'ld.  A.  (1993).  Efficiency  of  sampling  inspection  plans  for  attributes. 

Bull.  Int.  Statist.  Inst.  40,  6CI-697. 

11.  Hale!,  A.  (1965).  Bayesian  single  sampling  attribute  plans  for  discrete  prior 

istributions.  Mat.  Fys.  Skr.  Dan.  Vid,  Sclsk,  3,  No,  2,  08  pp,  Munksgaard, 
Copenhagen. 

12.  Hald,  A.  (1965)  Single  sampling  Inspection  plans  with  specified  cTcceptancc 
probability  and  minimum  costs,  Skand.  Aktuart Idokr.  4C.  (Forthcoming). 

13.  Half!,  .9.  'nd  Rou9)paard  F.  (19  3).  Some  limit  theoreme  for  the  Dodge-Romig 
/.(0:h  single  sampling  inspeccion  plans.  Sankliyn  ..  25,  255-60. 

Lehmann.  E.L.  (1958),  Significance  level  .and  power,  Ann,  Math.  Statist,  29, 
1167-1176. 


J 


-  33  - 

15,  Markbtfck.N,  (1950).  Bcll-Kontroll.  labeller  f*6r  cnkel  Provtognlng. 

Teknlskt  Meddelando  Nr.  7a  V.  Gvcrlgcs  Mekanf'orbund^  Stockholm. 

16,  Military  Standard  105  D  (IS63).  Sampling  procedures  an<l  tables  for  inspection 
by  attributes.  64  pp.  U.S.  Govern.  Printing  Office,  Washington  D.C. 

17,  Weibull,  I.  (1951).  A  method  of  determining  inspection  plans  on  an  economic 


basis.  Bull,  Int,  Statist.  Inst.  33,  OS-lCV*. 


